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We develop a new analytical method for solving real time evolution problems of quantum many- 
body systems. Our approach is a direct generalization of the well-known canonical perturbation 
theory for classical systems. Similar to canonical perturbation theory, secular terms are avoided in 
a systematic expansion and one obtains stable long-time behavior. These general ideas are illustrated 
by applying them to the spin-boson model and studying its non-equilibrium spin dynamics. 



The theoretical investigation of non-equilibrium quan- 
tum many-body systems has recently become a very ac- 
tive field of research due to seminal experiments in ul- 
tracold atomic gases (for example collapse and revival 
phenomena electronic nanostructures (for example 
transport beyond the linear response regime [2j|) and gen- 
erally qubit dynamics in the presence of quantum dissipa- 
tion. While non-equilibrium classical systems have been 
long studied, quantum systems in non-equilibrium hold 
the promise of many new phenomena yet to be discov- 
ered. On the theoretical side, progress is hindered by the 
notorious difficulty of solving non-equilibrium quantum 
many-body problems. Motivated by the recent experi- 
ments, significant progress has been made with powerful 
numerical methods like the time-dependent density ma- 
trix renormalization group (TD-DMRG) or the time- 
dependent numerical renormalization group (TD-NRG) 
[3, HI]. However, there are few reliable analytical meth- 
ods available, especially for non-perturbative problems (a 
noteable exeception is the real time RG method [6(). 

A key problem for analytical calculations is the appear- 
ance of secular terms in time t that grow with some power 
of t. Secular terms appear naturally if one attempts 
a direct perturbative expansion, e.g. in the Heisenberg 
equations of motion for the observables. Even if secular 
terms are multiplied by a small coupling constant, they 
inevitably invalidate perturbation theory for large times 
even for small coupling constants and make it impossible 
to draw conclusions about the long-time behavior. This 
problem is also very well-known from classical mechan- 
ics, dating back to studies of planetary motion in pre- 
vious centuries. In the context of analytical mechanics, 
its solution using canonical perturbation theory is well- 
established and can be found in any textbook (see, for 
example, Q)- The basic idea is to first transform the 
Hamiltonian to normal form using a canonical transfor- 
mation. One can then easily solve the equations of mo- 
tion for the new position and conjugate momentum vari- 
ables. Only after integrating these equations of motion 



does one reexpress the old variables in terms of the new 
time-evolved variables. It is well-established that this 
yields a much improved long-time solution without any 
secular terms even if the canonical transformation itself 
is only done perturbatively . Surprisingly, to the best of 
our knowledge to date no attempt has been made to im- 
plement an equivalent scheme based on unitary perturba- 
tion theory for quantum many-body systems. However, 
one key difference to classical systems is that in quantum 
many-body systems one is often dealing with a continu- 
ous energy spectrum, which makes naive unitary pertur- 
bation theory impossible due to vanishing energy denom- 
inators. A way to solve this specific problem has been es- 
tablished recently by means of the flow equation method 
@, Q (for related ideas see also the similarity renormal- 
ization scheme The central idea of the flow equa- 
tion method is to diagonalize a many-particle Hamilto- 
nian through a sequence of infinitesimal unitary transfor- 
mations that eliminate interaction matrix elements with 
large energy difference first before dealing with smaller 
energy differences. In this way one both reorganizes a 
perturbative expansion in an RG-like manner, which al- 
lows one to recover non-perturbative energy scales, and 
one avoids the above small energy denominator problem 
even for a continuous energy spectrum. 

In this Letter we develop the general framework for 
applying the flow equation method to analytically solve 
real time evolution problems in quantum many-body sys- 
tems in exact correspondence to canonical perturbation 
theory in classical mechanics. We will see that likewise 
secular terms are avoided and that one can obtain reliable 
results about the long-time dynamics even in a perturba- 
tive framework. We will then illustrate our approach by 
studying the real time evolution of the spin-boson model 
with an initially polarized spin and a relaxed bath. The 
spin-boson model is the paradigm of dissipative quantum 
systems and its non-equilibrium behavior has rece ntly 
been investigated using the TD-NRG method 111, |l2j, 
which motivates our choice. 
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Let us briefly review the basic ideas of the flow equa- 
tion approach (for more details see [§]). A many- 
body Hamiltonian H is diagonalized through a sequence 
of infinitesimal unitary transformations with an anti- 
hermitean generator fy(-B), 



B = 



dH(B) 
dB 



[f)(B),H(B)] 



(1) 



with H(B = 0) the initial Hamiltonian. The "canon- 
ical" generator [8| is the commutator of the diagonal 
part Hq with the interaction part Hi n t of the Hamilto- 
nian, i](B) = f [Hq(B), i/i„t (£?)]. Under rather general 
conditions the choice of the canonical generator leads 
to an increasingly energy-diagonal Hamiltonian H(B), 
where interaction matrix elements with energy trans- 
fer AE decay like exp(— B AE 2 ). For B = oo the Hamil- 
tonian will be energy-diagonal and we denote parameters 
and operators in this basis by ~ e.g. H = H(B = oo). 

The key problem of the flow equation approach is 
generically the generation of higher and higher order in- 
teraction terms in {T]), which makes it necessary to trun- 
cate the scheme in some order of a suitable systematic 
expansion parameter (usually the running coupling con- 
stant). Still, the infinitesimal nature of the approach 
makes it possible to deal with a continuum of energy 
scales and to describe non-perturbative effects. This had 
led to numerous applications of the flow equation method 
where one utilizes the fact that the Hilbert space is not 
truncated as opposed to conventional scaling methods. 
Examples are the evaluation of correlation functions on 
all energy scales in equilibrium problems [sj and non- 
equilibrium problems, where one cannot focus on low- 
energy degrees of freedom anyway (see, for example, the 
time-dependent Kondo model [l3| or the Kondo model 
with voltage bias [3). 

We will now utilize these features to develop an ana- 
logue of canonical perturbation theory in classical me- 
chanics for quantum many-body problems. The general 
setup is described by the diagram in Fig. [IJ where Y$>i) 
is some initial non-thermal state whose time evolution 
one is interested in. However, instead of following its full 
time evolution it is more convenient to study the real 
time evolution of a given observable A. This is done by 
transforming the observable into the diagonal basis in 
Fig. [T] [forward transformation) : 



dO{B) 
dB 



= [r,(B),0(B)} 



(2) 



with the initial condition 0(B = 0) = A. The cen- 
tral observation is that one can now solve the real time 
evolution with respect to the energy-diagonal H exactly, 
thereby avoiding any errors that grow proportional to 
time (i.e., secular terms): this yields A(t). Now since 
the initial quantum state is given in the B = ba- 
sis, one undoes the basis change by integrating ^ from 
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A{t) = e^Ae-* 6 ', |*,) 



FIG. 1: The forward-backward transformation scheme in- 
duces a non-perturbative solution of the Heisenberg equations 
of motion for an operator. U denotes the full unitary trans- 
formation that relates the B = to the B = oo basis. Ilal 



B = oo to B — (backward transformation) with the 
initial condition 0(B — oo) = A(t). One therefore effec- 
tively generates a new non-perturbative scheme for solv- 
ing the Heisenberg equations of motion for an operator, 



A(t) = e tHt A(0) 



-iHt 



in exact analogy to canonical 



perturbation theory. Notice that it is the last step of the 
backward transformation that distinguishes this scheme 
from the flow equation evaluation of equilibrium correla- 
tion functions [9(: The equilibrium ground state or ther- 
mal states are in fact more easily expressed in the B = oo 
basis (since H is energy-diagonal) than in the B = (in- 
teracting) basis. It should be mentioned that the same 
forward-backward transformation scheme with respect to 
some given initial quantum state has recently also been 
successfully employed by Cazalilla [TJ| for studying the 
nonequilibrium Luttinger model. The main difference to 
our approach is that the bosonized Luttinger Hamilto- 
nian becomes quadratic, which makes it possible to work 
out the unitary transformation exactly in (th e same 
holds in [l]|): therefore stability questions regarding sec- 
ular terms for a generic interacting system do not arise, 
which are the main focus of our work. 

We now illustrate the general idea of our approach 
by studying the spin-boson model, which serves as a 
paradigm in dissipative quantum physics and for qubit 
dynamics (for a review see, for example, [l3|) : 

H = -y a x + \a z £ X k (bl + b k ) + ]T u k b\b k . (3) 

k k 

It describes a two state system coupled to a bath of 
harmonic oscillators. The effect of this dissipative en- 
vironment is encoded in the spectral function J(w) *= 
J2k ^k ^(^ ~ w fe)- I n t ne sequel Afe is considered a small 
expansion parameter. In this Letter we will only study 
the zero temperature case, T = 0, although the gener- 
alization to nonzero temperature is straightforward. We 
use the following generator for the unitary flow [TJj]: 



V (B) = ia^V^ib. + bb + o-^rthh-bt) 

k k 

+ ■■(h + b{)(h-bl):, (4) 



k.l 



with independent coefficients: 
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Normal-ordering is denoted by : ... : , which serves as 
a systematic scheme to truncate the infinite sequence of 
higher and higher operators generated by (p}. Higher 
normal-ordered terms than the ones contained in ([3]) are 
neglected in the flow of the Hamiltonian, which amounts 
to neglecting small (of order X k ) higher order cumu- 
lants in the Hamiltonian (this approximation is reliable 
for any super-Ohmic bath and for an Ohmic bath with 
a < 0.2, for more details see @, [Hj]). If one is interested 
in equilibrium properties, normal-ordering is performed 
with respect to the equilibrium ground state, b k b\, =: 
b k b\, : +5 kk > n(k), where n(k) is the Bose-Einstein dis- 
tribution. However, later we will be interested in the 
real time evolution of a non-thermal initial state \^>i). 
Hence, in order to minimize our truncation error, we 
write more generally b k b\, =: b k b\, : +S kk > n(k) + Ckk', 

where C kk > =' | b k b\, | - 5 kk > n(k). The flow of 
H(B) generated by this 77 is 

dB V ^fc + A w 

^ = -K-A) 2 Afe + 2^ % ,A z . (7) 

The derivation of ([6]) and is discussed in detail in 0, 
[Til ]. The diagonalized Hamiltonian for B — 00 is 

H = -^a x + Y,"kblb k , (8) 

k 

where A = A(B = 00) is the renormalized tunnel- 
ing matrix element. For example for an Ohmic bath, 
J(u>) = 2aujQ(uj c — u>), the renormalized tunneling ma- 
trix element derived from the solution of the flow equa- 
tions 0, [3] has the correct non-perturbative behavior 
0, Aoc A (A/u; c ) a/1 - a . 

The observables in the B = 00 basis are given by solv- 
ing |(3J) for a suitable ansatz for the flowing observable 
Q. For example, 

a x (B) = h(B)a x +a z J2(xk(B)b k +Xk(B)bl^ 

k 

+ a(B) + ia y (MB) b k - Ji k {B) b[) (9) 
fc 

where higher normal-ordered terms generated in 0(X k ) 
during the flow ||3J| are again neglected. The differential 
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FIG. 2: Real time evolution of the spin expectation value 
(&z(t)) starting from a polarized spin in ^-direction with a 
relaxed Ohmic bath (see text) for two different values of a 
and Wc/A = 10. The full lines are the flow equation results, 
the dashed lines TD-NRG curves for A = 2.0 and the dotted 
lines for A = 1.41. The TD-NRG results are courtesy of F. 
Anders, see [llj. The various curves agree extremely well 
except for very long times shown in the insets. 

equations describing this flow take the following form: 

% = -J2{^ ) ^k + x k ) + Vk Z) ^k + fik)) 

k 

-4E# )C «(X/+*,) 

k,l 

= 2hr>M +J2( %l ( Xl + Xl)+V lk (Xi-Xi)) 
% = E H y) ^k + H) + vl z) (Xk + x k i) , (10) 

k 

with the initial conditions h{B = 0) = l,Xk{B = 0) = 
fi k (B = 0) = a(B = 0) = 0. For A 6 supp J(lu) the 

observable decays completely, h = h(B = 00) = 0, im- 
plying decoherence The ground state expectation 

value of o x is then given by a = f a(B = 00). 

For real time evolution problems we now solve the 
Heisenberg equations of motion in the diagonal basis, 
a z (t) — e lHt <r z e~ lHt . The result is straightforward 

Xk (t) = (xfe(O)cos(At) +i£ fc (0)sm(At))e- <w * t (11) 
Afc(*) = (/i fe (0)cos(At)+ixfc(0)sin(At))e" M '= t , 

while h and a remain unchanged. In complete analogy to 
canonical perturbation theory, we next undo the unitary 
transformation ||3J|. The values of h, a, Xk(t), p-kif) are 
used as initial values in the system of differential equa- 
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FIG. 3: Ground state expectation value of a x : Comparison of 
flow equation results (curves) and NRG data (squares) from 
for an Ohmic bath with damping a. The results are for 
u> c /A = 25, 28.6, 33.3, 40, 50, 66.7, 100 from top to bottom. 



h(t) a x +a z ^ (xk(t) b k + xk(t) b f K 

k 

a(t)+ia y J2Uk(t)b k -fl k (t)bi) (12) 



tions (10)) at B — oo, which is then integrated back- 
wards to B = 0. This yields h(t), a(t), Xk(t), Hk(t), which 
parametrize the time-evolved operator o~ x , 



cr x (t) 



in the original basis of the problem. Thereby the 
forward-backward transformation scheme induces a non- 
perturbative solution of the Heiscnberg equations of mo- 
tion, compare Fig. [T] 

For the purposes of this Letter, we focus on the nu- 
merical solution of the above differential equations by 
discretizing the bosonic bath with O(10 3 ) modes (notice 
that an approximate analytical treatment is equally pos- 
sible). The initial quantum state is taken as spin 
up, {^il&zl^i) — +1, with a relaxed bath with respect 
to this fixed spin. This yields Ckk' = AfcAfc/ /AujktOk' ■ We 
have implemented the numerical solution for all compo- 
nents of the spin degree of freedom. In order to assess the 
accuracy of our approach, the time evolution of {o~ z {i) is 
shown in Fig. [2] and compared with TD-NRG data for two 
values of the discretization parameter A. One finds excel- 
lent agreement except for very long time scales (shown 
in the insets of Fig. [2|), where the TD-NRG discretiza- 
tion error becomes noticeable (since the curves depend 
on A). The flow equation solution for the observable 
shows that it approaches its flow equation equi- 
librium expectation value (<j x )qs with an absolute error 
below 10 -2 for long times. A comparison of (c x )gs with 
exact numerical results using NRG [l9[ in Fig. [3] again 
shows very good agreement. 

Summing up, we have shown how to implement an 
analogous scheme to canonical perturbation theory for 
quantum many-body systems. Using a simple but non- 
trivial example, we could demonstrate that the well- 
established advantages of canonical perturbation theory 



versus naive perturbation theory carry over to our uni- 
tary perturbation approach as well, in particular the ab- 
sence of secular terms in real time evolution problems. 
Our results are stable in the long-time limit (see Figs. l2|3( l 
and can be improved systematically in a uniform manner 
(as a function of time) by higher orders of the calculation. 
The underlying scheme of infinitesimal unitary transfor- 
mations permits to study non-perturbative effects [20j ]. 
Similar to the role of canonical perturbation theory in 
analytical mechanics, our approach should be useful for 
other real time evolution problems from impurity systems 
to lattice models in quantum many-body physics 21 1 . 
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